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Abstract 

We calculate the low-temperature series of the free energy in absolute- 
value solid-on-solid (ASOS) model to order u^"^ using finite-lattice 
method. The property of the obtained series and the behavior of 
their Pade approximants confirms us that the roughening transition 
in ASOS model is of Kosterlitz-Thouless type. 
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1 Introduction 



Recently the low-temperature series in d = 3 Ising model or equivalently 
strong-coupling series in d = 3 Z2 lattice gauge theory have been extended 
to higher orders using finite-lattice method 0, |^. This method avoids the 
problem involved in the graphical method, in which it is rather difficult to 
give the algorithm for listing all the diagrams up completely that contribute 
to the relevant order of the series. The low-temperature series calculated by 
finite-lattice method in d = 3 Ising model are for the true inverse correlation 
length (which is equivalent to the mass gap in lattice gauge theory) free 
energy ||^, 0, magnetization and zero-field susceptibility 0, surface tension 
(which is equivalent to the string tension in lattice gauge theory) and 
the second moment of the correlation function f^. We should mention that, 
using a modification of the shadow-lattice technique, low-temperature series 
were obtained for free energy, magnetization and susceptibility which are 
longer than those of reference p, ^ . The second moment was also calculated 
by the same technique [0, in which exact series is shorter and estimated 



series is longer than that of reference 

As for the surface tension we 0] calculated the low-temperature series 
to order u^'^ with u, . = exp (—4:3. . ) and found that the series coeffi- 

laing Ising -t^ \ ' I sing/ 

cient change their sign at the order of u^^ . The surface tension in three- 

^ Ising 

dimensional Ising model or the string tension in three- or higher-dimensional 



lattice gauge theory suffers from the roughening transition [O, O, 13] and 



it is expected to exhibit Kosterlitz-Thouless type singularity like 



f{u) = A{u) exp [-c{ur - u)-^/^] + B{u). (1) 

It has the essential singularity at the roughening transition point Ur- As was 
pointed out by Hasenbusch and Pinn |T5[, the sign-change of the coefficients 
is just the signal of the K-T type singularity. In fact if we expand the function 
(|ip in terms of u we would obtain a series with sign-change. The order of 
the sign-change depends on c and Ur- The Pade approximants of the low- 
temperature series for the surface tension appear to converge below the order 
of the sign-change [|1^. They, however, start to scatter just at the order of the 



sign-change and appears to converge again for higher orders [0. The Pade 
approximants of the maximal-order series for the surface tension exhibit good 
convergence but their values overshoot the high precision Mote Carlo results 



[O. We think that our series are yet too short for the asymptotic regime to 
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be reached. Unfortunately it would be difficult to extend the series further by 
the presently available technique of the finite-lattice method and computer 
resources. 

Absolute-value solid-on-solid (ASOS) model is an approximation of the 
interface of d=3 Ising model. Overhangs and disconnected parts are ne- 
glected in the former. The ASOS model is also expected to exhibit K-T 
type phase transition. Its free energy corresponds to the surface tension of 
d = 3 Ising model and is expected to behave like (|l]). Hasenbusch and Pinn 
[ [Tsll calculated the low-temperature series for the free energy of this model 
to order -u^^ using the finite-lattice method, extending the previous series 
by Weeks et al [|l^, and they found the expected sign-change of the coeffi- 



cients at the order of u^^. They analyzed the series by Pade approximants, 
compared with their Monte Carlo data and obtained the results similar to 
the case of the surface tension in d = 3 Ising model. They concluded that 
the Pade approximants are inconsistent to the Monte Carlo data. Only two 
terms after the sign-change, however, would be too short to make a definite 
conclusion. 

In this paper we extend the low-temperature series for the free energy 
of ASOS model to order using the finite-lattice method. An improved 
technique is introduced to extend the series, which can also be applied to 
the low-temperature expansion of the spin systems whose spin-variable takes 
more than 2 values such as the spin-S* Ising model with S" > 1 or the g-state 
Potts model with q > 3. The obtained longer series enables us to make a 
definite conclusion on the roughening transition in the ASOS model and on 
the convergence of the Pade approximants for the free energy. 



2 Algorithm 

Here we give the improved algorithm of the low-temperature expansion for 
SOS model that enables us to obtain longer series than the case when the 
finite-lattice method was applied naively. Let us consider two-dimensional 
Lx X Ly rectangular lattice Aq. The free energy density in the infinite- volume 
limit is given by 

/ = - lim -i- In (2) 
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with Hamiltonian 

Z(L„ L,; = ^exp ( - /5 ^ \hi - hj\) (3) 

{ft} <ij> 

where the variable hi at each site i is restricted to —H^ hi < 
{Hj^^H^ > 0). The low-temperature series is calculated with respect to 
the expansion parameter u = exp(— 2/5). We take the boundary condition 
that all the variables outside Aq are fixed to be zero. 

Let us consider the set {A} of all two-dimensional rectangular sub-lattices 
of Aq. The sub-lattice A is denoted by its size x ly and its position in Aq. 
We define of A as 

H(A; h+, = - In [Z(A; h+, ]. (4) 

In the calculation of the partition function Z{A; h^,h_) the variable hi inside 
the A is restricted to the range of — /i_ < hi < h^ /i_ > 0) and all the 
variables outside are fixed to be zero. We define W oi A recursively as 

W{A; h+, h.) = H{A; h+, h.) - ^^'^ W{A'; h'^, h'_). (5) 

A'CA, 
h'j^<h+,h'_<h- 

Here the summation X^^'^ implies that the VF(A; h+,h_) should be excluded in 
the summation. We note that H{A; and 1^(A; is independent 

of the position of A. We know 

H{Ao-H+,H^) = W{A;h+,h^) 

ACAO, 
h+<H+, h-<H_ 

J2 {L,-h + l){Ly-ly + l) 

h+<H+, h-<H- 

xW{l:c X ly;h+,h-). (6) 

Taking the infinite-volume limit we obtain 

/ = lim -^H{Ao;H^,H_) 

= ^ W{l,Xly-h+,h^). (7) 
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In the standard cluster expansion of the free energy for this model, a 
cluster is composed of polymers and each of the polymers consists of sites 
that are connected by nearest-neighbor bonds. A value hi{^ 0) is put to each 
site i of the polymer. We can assign to each cluster two numbers hp and hm 
that are the maximum of hi{> 0) and —hi{> 0), respectively, in all the sites 
of the polymers the cluster consists of. Then, we can prove that the Taylor 
expansion of W{lx x ly;h+,h_) with respect to u includes the contribution 
from all the clusters of polymers in the standard cluster expansion that have 
hp = h+ and hm = h^ and that can be embedded into Ix x ly lattice but cannot 
be embedded into any of its rectangular sub-lattices. The series expansion of 
W{lx X ly] h+, hJ) starts from the order with n = lx + ly + 2(/i+ -|- hJ) — 2. 
So, to obtain the expansion series to order m^, we should take into account 
all the rectangular lattices that satisfy Ix + ly + '^{h+ + h_) — 2 < N in the 



u 



N 



summation of eq. (|^) and evaluate each of the W^s to order 

All the above algorithm can be coded on the computer and most of the 
CPU time and main memory are used for the evaluation of the relevant 
partition functions Z{lx x ly]h+,h-) to order . The partition functions 
are calculated using the transfer matrix formalism based on building up finite- 



size lattices one site at a time [18]. The necessary CPU time and memory 
are proportional to Nlxly {h+ + h_ + and + /i_ + 

respectively. In table |^ we list the maximum Mimp{N) of this factor 
A^ + /i_ + i^^^^C-^M for the necessary memory in calculating the series 
to order by this improved algorithm. 

It is quite natural to introduce the partition function for the finite-size 
lattice with the range of the variable hi restricted as in equation (^. It comes 
from the analogy with the algorithm of the low-temperature expansion for 
the surface tension in d = 3 Ising model by the finite- lattice method . We 
should stress that it enables us to calculate much longer series than the case 
when we would apply the finite-lattice method naively by considering the 
ASOS model as a two-dimensional spin system. In the latter case we would 
define 

if(A) = -ln[Z(A)] (8) 

and 

W{A) = H{A) - J2 W{A'). (9) 

A'cA 

for the rectangular x ly lattice A. Here the variable hi takes all integer 
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values — oo < hi < +00 in the calculation of the partition function Z{A). 
Then the free energy density in the infinite- volume limit would be given by 

f=Y^W{hxly). (10) 

For the expansion series to order m^, we should take into accout all the rectan- 
gular lattices that satisfy Ix + ly < N m the summation of equation (p!0|). The 
necessary CPU time and memory to evaluate each Z{K) to order for x ly 
lattice is proportional to M^./j^(2[A^/2] + l)™'^ and iV(2[A^/2] + l)°^''^ 
respectively. In table |l| we also list the maximum Mnaive{N) of this factor 
N{2[N/2] + 1)°^*° for the necessary memory in calculating the series to 
order by this naive algorithm. We can see that the necessary memory in 
the improved algorithm is much smaller than in the naive algorithm. 

We note that the improved algorithm above can be used in the finite- 
lattice method of the low-temperature expansion for the spin systems whose 
spin variable takes more than 2 values such as the spin-S* Ising model with 
S' > 1 or the g-state Potts model with g > 3. 



3 Series 

We have calculated the expansion series for the free energy / of ASOS model 
to order u"^^. The calculation was performed on FACOM-VP2600 at Kyoto 
University Data Processing Center. Used main memory was about 80Mbyte 
and total CPU time was about 25 minutes. The obtained coefficients for the 
free energy are listed in table 0, where a„ is defined as 

/ = E«n^"- (11) 

n 

We have checked that each of the W{lx x Z^,; /i_)'s in equation starts 
from the correct order in u discussed in the previous section. The series to 
order u^"^ for ASOS model coincides with those obtained by Hasenbusch and 



Pinn |T5| and we have obtained new 11 terms. Also listed are the coeffi- 
cients for the surface tension a in d = 3 simple cubic Ising model obtained 
previously by the author [0. The coefficients are defined by 



cr 



2/3/..„, + E (12) 
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where uismg = exp (— 4/?/si„g) is the low-temperature expansion parameter in 
the Ising model. 



4 Series Analysis 

It is to be noticed that the coefficients of the free energy in ASOS model 
are all negative to order u^^ and their sign changes at the order of u^^ [p!5|] . 
We find no another change of sign in higher orders to u'^^. This property of 
the sign-change is common to the surface tension in ci = 3 Ising model. As 
was pointed out in reference |]15| it supports the claim that the roughening 
transition of these models is of Kosterlitz-Thouless type. Let us consider a 
K-T type function and its expansion in u as 

/^^(u) = Au^ exp [-c{ur - = ^ ^u" (13) 



and fit this expansion to the expansion of the free energy ([IT]) of ASOS 
model, so that the ratios of their successive coefficients coincide as an/cin-i = 
a^^ /a^'^^ for n = 11 and n = 23, respectively. The order = 11 is just 
where the change of sign occurs and n = 23 is the highest order of our 
series. Then we obtain Ur = 0.214 and c = 0.533. In table ^ we list the 
ratios a„/a„_i for the free energy and the fitted K-T type function. We 
note that the ratios agree with each other within 3 per cent for the range of 
n = 6 to 23. Furthermore the ffited value of Ur is not far from the predicted 
roughening transition point for ASOS model as Ur = 0.207(9) from the series 
analysis of the surface width |l9l and Ur = 0.1994(1) from the Monte Carlo 



renormalization-group analysis |]20|. The difference between the ratios for 
n < 5 would come from the regular background term in the free energy, 
whose contribution can be expected to become negligible rapidly for higher 
orders compared with the contribution from the singular term. 

Next we will try the Fade approximation for the free energy f{u) and 
the internal energy E{u) = df/d(3 = —2uf'{u), although the free energy 
is expected to behave like equation ([l|) and one might think that the Fade 
approximation would not be appropriate for such a function. The authors 
of reference stressed the discrepancy between the Fade approximants of 



the internal energy to order u^"^ and their Monte Carlo data. Their series 
are, however, too short to give a definite conclusion. 
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In figures 1 and 2 we show the truncated series of the free energy and its 
Pade approximants for each order n of the series m. (3 = 0.85(m = 0.1827) and 
13 = 0.81 (m = 0.1979), respectively. The latter point of /? is very close to the 
predicted roughening transition point ||2^ . As for the Pade approximants we 
plot the [L, M] approximants {L+M = n-2) with (ri-2)/4 < L< 3(n-2)/4. 
Defective approximants that have poles on the real axis of u in the range of 
< u < 0.205 are excluded. The Pade approximants appear to converge 
once to some value for the shorter series. They begin to scatter just around 
the order of the sign-change and they converge again for higher order series 
but to another value. The convergence is so good in higher orders. The 
convergence of the truncated series is, on the other hand, much slower than 
the Pade approximants. In figure 3 we also give the behaviour of the Pade 
approximants together with the truncated series of the free energy in f3 for 
the highest order series. From these figures we can conclude that the Pade 
approximation works well for the long enough series of the free energy of 
ASOS model even very close to the predicted roughening transition point. 
This situation of the Pade approximants for the free energy of ASOS model 
is quite similar to the case for the surface tension in c? = 3 Ising model [l^, 0], 
in which presently available series has only 5 terms after the sign-change. We 
can expect that the higher order series for the surface tension in d = 3 Ising 
model will also give Pade approximants converging to the Monte Carlo data 
and the present discrepancy [0 will be resolved. 

In figures 4 and 5 we show the truncated series of the internal energy 
and its Pade approximants for each order n of the series for f3 = 0.85 and 
P = 0.81, respectively. We see that the behaviors of the truncated series and 
its Pade approximants in the order n are similar to those for the free energy 
given in figures 1 and 2. The convergence of the Pade approximants is good 
for higher orders, although the convergence for (3 = 0.81 is slower than in the 
case of /? = 0.85. On the contrary, the convergence of the truncated series 
is much slower than the Pade approximants. In figure 6 and table ^ we also 
give the behaviour of the Pade approximants for the internal energy versus 
f3 for the highest order series. The Pade approximants for the highest order 
series are consistent with the Monte Carlo data |15| even at /3 = 0.81 very 
close to the roughening transition point. 

One's question might be why the naive Pade approximation, which is the 
ratio of two polynomials, is so good even close to the roughening transition 
point in spite of that the free energy is expected to have essential singularity 
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as in equation (|T]). To clarify the reason we examine the Pade approximation 
for a test function 

g{x) = exp [-(1 - x)"-^/^] + exp (x). (14) 

The first term has the essential singularity at x = 1 and the second term is a 
regular background. In figures 7 and 8 we plot the exact function, the 23rd 
order truncated series and its Pade approximants of the test function g{x) 
and of its derivative xg (x), respectively. As for the test function itself, we 
can see that the Pade approximants exhibit so good convergence and give 
the values that are consistent with the exact function for x < 0.98. As for 
the derivative, we observe that the Pade approximants converge and follow 
the exact function well for x < 0.95. The derivative shows worse convergence 
than the test function itself, which would be because the derivative has the 
power singularity multiplied by the essential singularity. The truncated se- 
ries, on the contrary, starts to deviate from the exact function much earlier 
than the Pade approximants. We also plot in figures 9 and 10 the exact 
value, truncated series and its Pade approximants of the free energy and the 
internal energy, respectively, for each order n of the series in x = 0.98. The 
situation of the test function in figures 7 and 9 is quite similar to that of 
the free energy of ASOS model in figures 3 and 2 and the situation of the 
derivative in figures 8 and 10 is quite similar to the internal energy of ASOS 
model in figures 6 and 5, respectively. We have examined several other test 
functions that have the same structure as equation (|l]) and the situations 
are all similar to the case of the test function g{x). We can conclude that 
the Pade approximants of the function like equation (|I|), which involves an 
essential singularity, converges to the correct value even very close to the 
singularity point, if we use the expansion series in high enough order. 

5 Summary 

We have calculated the low-temperature series for the free energy of ASOS 
model to order using the finite-lattice method. An improved technique 
has been used to extend the series, which can also be applied to the low- 
temperature expansion of the spin systems whose spin-variable takes more 
than 2 values. The property of the obtained series and the behavior of their 
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Pade approximants have confirmed us that the roughening transition in the 
ASOS model is of Kosterlitz-Thouless type. 
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Table 1: Comparison of the factors Mimp{N) for the improved algorithm 
and Mnaive{N) for the naive algorithm which the necessary memory is pro- 
portional to in calculating the series to order . 



N 




Mnaive{N) 


2 


4 


6 


3 


6 


9 


4 


16 


100 


5 


20 


125 


6 


54 


2058 


7 


63 


2401 


8 


216 


52488 


9 


243 


59049 


10 


810 


1610510 


11 


891 


1771561 


12 


3072 


57921708 


13 


3328 


62748517 


14 


14336 


2392031250 


15 


15360 


2562890625 


16 


65536 


111612119056 


17 


69632 


118587876497 


18 


294912 


5808378560022 


19 


311296 


6131066257801 


20 


1562500 


333597619564020 


21 


1640625 


350277500542221 


22 


8593750 


20961814674106394 


23 


8984375 


21914624432020321 


24 


46875000 


1430511474609375000 


25 


48828125 


1490116119384765625 
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Table 2: The low-temperature expansion coefficients a„ for the free energy 
of ASOS model and 6„ for the surface tension of 3-d Ising model. 



Tl 


0"n 


On 


9 


n 
— Z 


n 
— Z 


o 
O 


A 

— 4 


— z 


A 


—in 


—in 





0/1 






— 1Q4/3 


—242/3 


7 


-172 


-150 


8 


-452 


-734 


9 


-3184/3 


-4334/3 


10 


-8862/5 


-32122/5 


11 


+1712 


-10224 


12 


+118604/3 


-106348/3 


13 


+284232 


+53076 


14 


+11450234/7 


+3491304/7 


15 


+42775116/5 


+74013814/15 


16 


+42439230 


+27330236 


17 


+203605872 


+160071418 


18 


+8589601858/9 




19 


+4397057716 




20 


+19982814764 




21 


+1885469273984/21 




22 


+4393383644538/11 




23 


+1760491870744 
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Table 3: Comparison of the ratio of a„/a„_i for the free energy and a^'^ / a^^^ 
for the fitted K-T type function f^^[u). The parameter for the K-T type 
function is Ur = 0.214 and c = 0.533. 



n 


'^n/'^n— 1 


"n /"n-1 


Q 
O 


Z.UUU 


— z.oyz 


A 

4 


z.oUU 


z.ioo 


K 

o 






6 


2.694 


2.658 


7 


2.660 


2.683 


8 


2.628 


2.591 


9 


2.348 


2.321 


10 


1.670 


1.634 


11 


-0.966 


-0.966 


12 


23.093 


22.345 


13 


7.189 


7.029 


14 


5.755 


5.635 


15 


5.230 


5.126 


16 


4.961 


4.869 


17 


4.798 


4.717 


18 


4.687 


4.618 


19 


4.607 


4.549 


20 


4.545 


4.499 


21 


4.493 


4.461 


22 


4.448 


4.432 


23 


4.408 


4.408 
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Table 4: Comparison of the Pade approximants of the low-temperature 
series and the Monte Carlo data in reference [15] for internal energy. The 
value of the truncated series is also listed. 



p 


0.85 


0.84 


0.83 


0.82 


0.81 


u 


0.1827 


0.1864 


0.1901 


0.1940 


0.1979 


truncated 


0.5917 


0.6252 


0.6590 


0.6919 


0.7222 


Pade [6/15] 


0.5886 


0.6202 


0.6509 


0.6797 


0.7060 


[7/14] 


0.5886 


0.6201 


0.6508 


0.6792 


0.7045 


[8/13] 


0.5886 


0.6202 


0.6509 


0.6795 


0.7054 


[9/12] 


0.5887 


0.6203 


0.6511 


0.6803 


0.7081 


[10/11] 


0.5887 


0.6203 


0.6511 


0.6803 


0.7079 


[11/10] 


0.5887 


0.6203 


0.6512 


0.6806 


0.7089 


[12/9] 


0.5887 


0.6203 


0.6513 


0.6807 


0.7094 


[13/8] 


0.5886 


0.6202 


0.6509 


0.6797 


0.7059 


[14/7] 


0.5886 


0.6202 


0.6509 


0.6796 


0.7059 


[15/6] 


0.5886 


0.6201 


0.6507 


0.6790 


0.7039 


MC(L = 256) 


0.5887(2) 


0.6205(2) 


0.6518(2) 


0.6811(2) 


0.7096(2) 


MC(L = 128) 


0.5884(4) 


0.6203(4) 


0.6516(4) 


0.6805(4) 


0.7095(3) 


MC(L = 64) 


0.5874(8) 


0.6498(7) 


0.6498(7) 


0.6795(6) 


0.7074(6) 
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Figure captions 



Fig. 1 : The values of the truncated series (+) and its Fade approximants 
(O) for the free energy versus the order n of the series in = 0.85(m = 
0.1827). The horizontal line represents the average of the Fade approximants 
for the highest order series of = 23. 

Fig. 2 : The values of the truncated series (+) and its Fade approximants 
(O) for the free energy versus the order n of the series in j3 = 0.81(m = 
0.1979). The horizontal line represents the average of the Fade approximants 
for the highest order series of = 23. 

Fig. 3 : The Fade approximants (O) together with the truncated series 
(dashed line) of the free energy versus u for the highest order series of n = 23. 

Fig. 4 : The values of the truncated series (+) and its Fade approximants 
(O) for the internal energy versus the order n of the series in f3 = 0.85. 
The horizontal line represents the result of the Monte Carlo simulation in 



reference [0] 



Fig. 5 : The values of the truncated series (+) and its Fade approximants 
(O) for the internal energy versus the order n of the series in f3 = 0.81. 
The horizontal line represents the result of the Monte Carlo simulation in 



reference |I5| ]. 



Fig. 6 : The Fade approximants (solid lines) together with the truncated 
series (dashed line) of the internal energy versus u for the highest order series 
of n = 23. The result of the Monte Carlo simulation (^)in reference is 
also plotted. 

Fig. 7 : The exact function (solid line), the 23rd order truncated series 
(dashed line) and its Fade approximants (O) of the test function g{x). 



Fig. 8 : The exact function (solid line), the 23rd order truncated series 
(dashed line) and its Fade approximants (O) of the derivative xg (x) of the 



test function. 



Fig. 9 : The values of the truncated series (+) and its Fade approximants 
{O) for the test function g{x) versus the order n of the series in x — 0.98. 
The horizontal line represents the exact value. 

Fig. 10 : The values of the truncated series (+) and its Fade approximants 
(O) for the derivative xg'{x) of the test function versus the order n of the 
series in x — 0.98. The horizontal hne represents the exact value. 
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